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Abstract. After a brief summary of the Newton-Cartan theory in a form which emphasizes its close analogy to 
general relativity, we illustrate the theory with selective applications in cosmology. The geometrical formulation 
of this nonrelativistic theory of gravity, pioneered by Cartan and further developed by various workers, leads to 
a conceptually sound basis of Newtonian cosmology. In our discussion of homogeneous models and cosmological 
perturbation theory, we stress the close relationship with their general relativistic treatments. Spatially compact flat 
models also fit into this framework. 



1 Introduction 



We hope that Klaus and Walter will accept this modest note as a tribute to their outstanding role as teachers 
of theoretical physics. In their courses they present not only elegant techniques and formal developements, 
but always emphasize the importance of basic concepts. "Rechnen kann jeder" , as Heitler used to say. In 
this spirit, we devote this article to a theme which is mainly of conceptual nature, and - as we hope - also 
of some pedagogical interest. 

We shall try to make it apparent that Newtons theory of gravity is much closer to general relativity 
(GR) than commonly appreciated. This has often been stressed in private conversations and letters by our 
inspiring teacher and colleague Markus Fierz. Here an example from a letter (Nov. 22, 1993): 



,,Es war um 1953, als ich meinen Newton-Aufsatz schrieb, dass ich Pauli sagte, auch in der allge- 
meinen Relativitatstheorie seien Raum-Zeit ,absolut', wie bei Newton. Darauf antwortete Pauli 
zu meinem Staunen: ,Sie verraten damit, das Grundprinzip der allgemeinen Relativitatstheorie: 
dass namlich Raum-Zcit-Matcrie nicht unabhangig voneinander gedacht werden konnen'. Wir 
konntcn uns dann aber einigen, indem ich zugab, dass hier das ,sine mutua actione' Newtons nicht 

'Dedicated to the sixtieth birthdays of Klaus Hcpp and Walter Hunziker (to appear in Helv. Phys. Acta) 
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gilt, obwohl im Ganzen diese Wechselwirkung klein ist. Was ich im Sinne hatte ist dies: Leibniz 
erklarte, der Raum sei nichts Wirkliches, sondern entspringe der Ordnung der Monaden (Kraft 
prastabilisierter Harmonie). Newton erwiderte hierauf: der Raum sei mehr als eine blosse Ord- 
nung. Denn der Abstand zweier Punkte habe einen Sinn ganz unabhangig davon, ob der Raum 
von etwas erfiillt oder leer sei. Newton legte also grosses Gewicht auf den metrischen Charakter 
des Raumes: dieser macht ihn zum Gegenstand der Physik, zu etwas Wirklichem ... " 

Every honest teacher of theoretical physics is confronted at a very early stage of a classical mechanics course 
with the following difficulty: After having introduced - in the spirit of L. Lange - the operational definition of 
an inertial frame, the question arises how to proceed when gravitational fields are present. In the traditional 
presentation of Newtons theory one maintains the fiction of an integrable (flat) affine connection, and puts 
gravity on the side of the forces, described by vector fields. A much more satisfactory formulation was given 
by Cartan and Friedrichs M . This denies the separate existence of a flat affine connection of space-time 
and a vector field describing gravitation, but puts gravity on the side of a more general dynamical connection 
which represents both inertia and gravitation. 

Historically, this important step of course was made first by Einstein when he created his general theory 
of relativity, but it is clearly independent of the relativization of time. Following Cartan and Friedrichs, 
numerous authors have elaborated on this idea. Here we mention only a selective list of contributions by 
Havas @, Trautman §, Ehlers §, |, |] and Kiinzle g §. 

In the first part of the present paper we give a brief summary of the Newton-Cartan theory, following 
mainly the work of H.P. Kiinzle, a former diploma student of Fierz at the ETH. Kiinzle's presentation, 
which uses the language of fibre bundles, appears to us as the most natural one, because it just replaces the 
role of the Lorentz group in GR by the Galilei group. Following this route, one arrives at a theory which is 
slightly more general than Newton's theory. The latter is only obtained after imposing a somewhat strange 
looking nonlinear condition for the Riemann tensor. The structural analogy of GR and the Newton-Cartan 
theory is, however, striking. In particular, the field equations look identical. 

In later sections we shall illustrate this also in more concrete terms with some selective applications in 
cosmology (homogeneous cosmological models and cosmological perturbation theory). This is perhaps not 
only an academic exercise, because much of the activity in cosmology, especially in connection with large 
scale structure formation, relies on the Newtonian approximation. We take this as a motivation for putting 
Newtonian cosmology on a conceptually firm basis. This has to be regarded as an extension of classic works 
by Heckmann and Heckmann and Schiicking jO]. One advantage which results is the possibility to 
choose the spatial sections as flat tori and thus describe compact cosmologies. 

All this confirms in more technical terms the remarks by Markus Fierz, quoted earlier. Einstein was wrong 
when he believed that his theory of gravitation incorporated the principle of Mach which is entirely in the 
spirit of Leibniz. This became already quite clear with the famous solution of Godel, but some relativists, 
notably Einstein himself, maintained the belief that Mach's principle might have something to do with the 
finiteness of space ^3|. That this is not the case was once and for all demonstrated by the "finite rotating 
universe" solution found by Ozvath and Schiicking . Space-time has really an independent existence and 
we are in fact still much closer to Newton than to Leibniz. 

A more detailed account of the material treated in this paper can be found in the diploma work by one 
of us JT|. 

2 Galilei spacetimes and their connections 

In what follows, M will always denote the space-time manifold and L(M) the principle bundle of linear frames 
with the structure group GL(4,IR). In GR space-time is endowed with a Lorentz metric g which defines 
a bundle reduction of L{M) to the orthonormal frame bundle 0{M) with the homogeneous Lorentz group 
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as the structure group. Conversely, each reduction of the structure group GL(4,IR) to the homogeneous 
Lorentz group gives rise to a Lorentz metric, because any element u £ L(M) over x S M can be regarded as 
a linear isomorphism of 1R 4 onto T X M , which maps the standard basis {e M } of IR 4 to the linear frame u. 

In a "nonrelativistic" gravity theory M has to be endowed with a Galilei metric, which consists of 
a one-form r and a symmetric semi-definit contravariant tensor field h of rank 3, satisfying h(-,r) = 
(Wtv — 0). The pair (7i, r) defines again a bundle reduction of L(M), this time with the homogeneous 
Galilei group as structure group. The reduced bundle consists of all frames {e v } in L(M), satisfaying 

r(e ) = l, h(0»,6°) = O, h{6\^)=5^ (i,j= 1,2,3), (2.1) 

where {8^} denotes the dual frames, (9^,e v ) = 6%. Since h{0^,T) = 0, these equations imply t = 9°. 

Conversely, a reduction of the structure group GL(A, IR) to the homogeneous Galilei group gives rise to a 
Galilei metric (h, r). This just reflects the fact that the homogeneous Galilei group (without time reflections) 
is the subgroup of GL(4,H) which leaves the standard Galilei metric of H 4 invariant. The latter is defined 
by equations ( |2.1| ) for the standard basis {e M } of 1R 4 and its dual. This defines the flat Galilei spacetime. 
With this notion it is also clear what is ment by locally flat Galilei spacetimes. These can be characterized 
as follows. 

Proposition 1 A GalileA spacetime (M,h,r) is locally flat iff the following two conditions are satisfied: 

(i) dr = 0; 

(ii) the induced Riemannian metrics on the integral manifolds defined by r are locally flat. 



From now on we shall only consider bundle reductions to the identity components G\ of the homogeneous 
Galilei group (orthochronous Galilei group), which we shall denote by Q (M, G+). The corresponding frames 
are then space and time oriented. 

It is now clear, how to define a Galilei connection on (M, h, r). This is a connection in the corresponding 
principle bundle G(M, G+), which we describe by a connection form u>, satisfying the usual conditions. There 
is a natural characterization of Galilei connections: 

Proposition 2 A linear connection T on a Galilei manifold (Af, h, r) is a Galilei connection iff 

V/i = 0, Vt = 0, (2.2) 

where V denotes the covariant derivative with respect to T. 



We consider only symmetric connections. For these the second equation in (2.2) implies dr — 0. Thus the 
distribution defined by the 1-form r is integrable. The corresponding maximal integral manifolds are the 
spatial sections of constant time. Vectors tangent to these sections are annihilated by r and are called 
spacelike (or horizontal). Tangent vectors which are not annihilated by r are called timelike. If t(V) = 1 
we say that V is a timelike unit vector. 

In contrast to Lorentz manifolds there is no unique symmetric Galilei connection on a Galilei manifold. 
It is instructive to see this in the light of a famous theorem by Weyl [[u] and Cartan jl5| . Since this is not 
so well-known (even among relativists) we state it here: 



Theorem 1 (Weyl, Cartan) For a closed subgroup G of GL(n,lR), n > 3, the following two conditions are 
equivalent: 



(i) G consists of all elements of GL(n, IR) which preserve a certain non- degenerate quadratic form of any 
signature; 
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(ii) For every n- dimensional manifold M and for every reduced subbundle P of L(M) with group G, there 
exists a unique torsion-free connection in P. 



It turns out that the set of symmetric Galilei connections is in 1:1 correspondence with the set A 2 (M) 
of 2-forms on M. The two equations (2.2) imply that the difference of two connection forms (Christoffcl 
symbols) is given by a tensor field of the following type 

Sap = 2T {a K mx h^, (2.3) 

where K Q( g are the components of a 2-form k. Special symmetric Galilei connections can be described as 
follows. Choose a timelike unit vector field V and define the covariant metric h b (relative to V) such that 
its components h^ LV satisfy 

K„v v = o, h^h x » = Tfi v u , (2.4) 

then 

V Kb = h^(V%h a(a T m + \h x %h Xa K ) - h p{a h^f3 - r (a V% (2.5) 

defines a symmetric Galilei connection. This is actually the unique symmetric Galilei connection which 
satisfies also 

V P V% = 0, h pa V% - h pB V a p = 0. (2.6) 



With respect to (2.5) the vector field V is then geodesic and rotation free. 



Note that relative to a Galilei frame with V = eo, equations (2.4) reduce to h^o — and h%j = 5ij\ 
thus h b = 5 lj e i i»9 j . This is a Riemannian metric on the leaves of the foliation defined by r. Clearly, the 
restriction of h b on an integral manifold is independent of V, because this is just the inverse of the restriction 
of the metric h. 

One can show that the integral manifolds (sections of constant time) are totally geodesic for any symmetric 
Galilei connection and that the induced connection on a leave coincides with the Levi-Civita connection 
corresponding to h b . 



The Newton-Cartan theory of gravity involves special symmetric Galilei connections of the form 

■a/3 — 1 a/3 ' °a/3 

where S£p is given by (2.3) with dn = 0. Such connections will be called Newtonian. 



Kb= V Kb + S^ (2.7) 



We need also a characterization of locally flat Galilei spacetimes. 



Proposition 3 For a Galilei manifold (M, h, r) with symmetric Galilei connection T the following state- 
ments are equivalent: 

(i) the Galilei manifold is locally flat; 

(ii) W v := h^Ph va R pa = 0; 

(iii) = ott^Tv) for some 1-form a. 

Recalling that any Newtonian connection can be expressed in terms of the Galilei metric (h,r), a timelike 
unit vector field V and a closed 2-form n, the question arises, when - for a given Galilei metric - a change 
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of k and V does not affect the Newtonian connection. One can show that this Newtonian gauge group 
is given by 



, (2-8) 



where / is a smooth function and w, A are 1-forms with k = jdA. (For an elegant proof see fl6|.) 

For many purposes it is useful to work in adapted coordinates: As a consequence of the Frobenius 
theorem for the integrable distribution defined by r, we can introduce local coordinates (t, x 1 , x 2 , x 3 ) in the 
neighborhood of any spacetime point such that r = dt and r(9j) = 0. The integral manifolds are then the 
slices of constant t (absolute time). Furthermore, the condition h(-, r) = implies that h = h^di ® dj. In 
adapted coordinates with (x° = t) 

r = dx°, h = h ll d l ®d ] (2.9) 

and the timelike unit vector field 

V = d , (2.10) 
the expressions for the Christoffel symbols of a symmetric Galilei connection become 



= o, rg b = h ac ( Kbc + i/i 6c , ) 

Tqq = 2h ab K b. V£ c = \h ad {hdb,c + hdc.b — hb c ,d)- 



(2.11) 



Here (h%j) is the inverse matrix of (h lJ ), in other words h — h^dx % ® dx 3 . The last equation in (2.11) 
proves our previous statement, that the induced connection on the slices of constant time is the Levi-Civita 
connection for the restrictions of h b . 

In addition to the space metric h we introduce the time metric g = t ® r. Clearly, 

gaphfr = 0. (2.12) 

In contrast to GR, the two (degenerate) metrics h and g are not the inverses of each other. 

3 The Newton-Cartan theory 

After these geometrical preparations we can now formulate the Newton-Cartan theory in a form which 
emphasizes its close analogy with GR. The theory consists of three parts: 

I Spacetime is a Galilei manifold (Af, h, r), with a Newtonian connection T. 

II Matter is described in part by a symmetric contravariant energy momentum tensor T a @ with vanishing 
covariant divergence (relative to T): 

V p T af} = 0. (3.1) 

III The field equations are 

R a0 = 8irG(T a0 - iffa/jT) - Ag a p, (3.2) 
where T a(3 := g acr g Si pT ap , T := g ap T ap . 
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In this formulation we have basically replaced the Lorentz group by the Galilei group. Several remarks are 
in order. 



First, it has to be emphasized, that (3.1) is not a consequence of the field equations. This is related to 
the fact that a Galilei metric does not fix the connection. 

The specialization to a Newtonian connection lookes somewhat mysterious. There is an equivalent for- 
mulation of this in terms of a symmetry of the Riemann tensor : 

a(R((3 t ,X)Y) = P(R(a t ,Y)X) (3.3) 

for any covectors a, f3 and vectors X,Y; (j denotes the map a cr = h(-,a). In index notation (3.3) reads 

K">R<$ pS = h^R} pp . (3.4) 

In GR, where h a P is the inverse of g a p, this symmetry is automatically satisfied. Since the Galilei metric 
does not fix the connection, we have the freedom to impose (|3.4|) as a further restriction. 



The field equations, which can also be written in the form (g a p = T a Tp) 

Raf3 = ^TlGpTaTp - k-T a Tp, p := T = T a TpT a(3 , 



(3.5) 



allow us to introduce Galilei coordinates: Clearly (3J5) implies R a/3 — and thus by Proposition || the 
Galilei manifold is locally flat. We can therefore specialize the adapted coordinate conditions ( |2.9[ ) even 
further such that 



dx v 



s ij di <g> di 



(3.6) 



In adapted coordinates we have Rij = as a consequence of the field equations, which also implies that the 
threedimensional time slices are locally flat. 



In Galilei coordinates the Christoffel symbols ( |2.11 ) simplify to 

o, r a 00 = 2h ac K 0c , r a Qb = h ac K bc , r a bc = o. (3.7) 

The Newton-Cartan theory is slightly more general than Newtons theory of gravitation. This can be seen 
by writing the field equations (3.5) for A = in Galilei coordinates. Inserting (3.7) one finds 



1 a/3 



2Kq j Kij K ' 



AirGp, 



(3.8) 



and 



s/j = 0. (3.9) 

In addition to this we also have dn = 0. We would obtain Newton's theory if the Galilei coordinates could 
be choosen such that Kij — 0. (Note that we can still perform time dependent rotations and translations.) 
Now, one can show Q that this is possible if and only if the following nonlinear condition for the Riemann 
tensor is imposed 



(3.10) 



h^"R^X p x = 0- 

Relative to Galilei coordinates which satisfy also Ky = 0, we obtain for g = — 2(/coi, K02, ^03) from ( |3.8| ) and 
dn = the basic equations of the Newtonian theory: 



div<7 = —AirGp, 



curly = 0. 



(3.11) 



Ehlers has shown g, that the strange condition (3.10) can be deduced from a spatial boundary condition 
at infinity which can naturally be imposed for the description of isolated systems. 



G 



One advantage of the geometrical formulation of the Newton-Cartan theory is that the spatial sections 
can also be chosen as flat tori. This enables us to describe spatially compact cosmological models. Some 
cosmological aspects will be presented later. 



Finally note that equation (3.8) reads (including the cosmological term) 

divg = -4nGp + A + k, 3 k ,j . (3.12) 
This shows that KijK 1 ^ acts (like a positive A) as a repulsive source. 



4 Fluid models in the Newton-Cartan theory 

This section serves mainly as a preparation for our later discussion of Newtonian cosmology. 

We introduce again a distinguished timelike unit vector field V on the Galilei manifold (M, h, r) with 
time metric g = t®t. The integral curves of V define a family of fundamental observers. Note that r(V) = 1 
translates into r a — g a pV^ ■ The matter model is assumed to be an ideal fluid with four velocity u, which is 
also a timelike unit vector field. We begin with some kinematical considerations which are familiar in GR. 

It is useful to introduce the projection operator P : T X M — > S x := kerr^ from the tangent spaces onto 
the horizontal (i.e., spacelike) subspaces definied by 

P(X) = X-g(X,V)V. (4.1) 

Clearly, P{V) — and r (P(X)) = 0. The components of P are 

PZ = %-9»\V X V V . (4.2) 

As before, ha U denotes the components of . We have the identities 



//,/• 

PY.=h^ x h x \ P£h x » = h^, P^h Xu = h^, 



pA r>v pf per pp o 



(4.3) 



For the covariant derivatives of u and V one verifies readily the following facts: 

VxV and Vxu are horizontal, (4-4) 
The vorticity (relative to V) is the skew symmetric bilinear form 



Q(X, Y) = \[h{V P{Y) u, P{X)) - h(\7 P{x) u, P(Y))} (4.5) 



and the (rate of) strain is 

Q(X, Y) = l -[h{V P(Y)Ul P{X)) + h(V P{X) u, P(Y))}. (4.6) 

The expansion rate is 

= h al3 e a0 (4.7) 
and the (rate of) shear is the trace-free part of the strain 

a(X, Y) - Q(X, Y) - ±6h(X, Y). (4.8) 
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While these quantities have the usual interpretation for the fluid motion relative to V, they are, unfortunatly, 
not tensoriel. They are, however, simply related to the contravariant tensor fields 17", with components 



(4.9) 



e 



Indeed, the components of fli.q ) and (4J3) are given by 



— h ap hp a il pa 



©a/3 — h ap hp a Q f 



With (O) and (;h4) one finds that 9 is simply given by 



and the covariant derivative of the velocity field can be decomposed as follows 

haXU^ = Q a/3 + fl a f3 + h ap V X U P x gf3 a V' 7 



(4.10) 



(4.11) 



(4.12) 



(4.13) 



With the help of (4.13) we can now derive a Raychaudhuri equation in the Newton-Cartan theory. As in 
GR we start from the identity 



which gives 



„.a a Tja a 



(4.14) 



With the help of (4.13) and the identities collected in ( |4.3| ) and (4.4) one can write the second term on the 
right as follows 



(4.15) 



The first term on the right in (4.14) is 

(u fi u? p ). a = div(V„u). (4.16) 
After a few steps (see Jl^]), we arrive at the following two equivalent forms of the Raychaudhuri equation 



div(V u u) = V„6» + -6 2 + h ap h pa {a pa a a[} - VL pa tt af} ) + Ric(u, u) 

= v„u + -e 2 + h ap h Pr7 (<T pa a a(3 - n pa n af3 ) + Ric(u, «). 

o 

Note that these equations hold for any Galilei manifold with a symmetric Galilei connection. 
At this point we use the field equations ( |3.5|) and obtain (with T a Tp — gaaQfipU^u 9 ) 

div(V„w) = V u 9 + V + h ap h^{a p<J a af} - fl p M a p) + <inGp - A. 

This equation will play an important role. 

Now we consider an ideal fluid with the energy momentum tensor 

T = pu® u + ph. 
From VT = one obtains the continuity equation 

div(pu) = 



(4.17) 



(4.18) 



(4.19) 
(4.20) 
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and the Euler equation 



-div (ph). 



(4.21) 



In contrast to GR, equation (4.2C) is a conservation law, because it is for a symmetric Galilei connection 
equivalent to (L u denotes the Lie derivative with respect to u) 



L u (pvol) = 0, 



(4.22) 



where vol is the standard volume r Avohj, V0I3 being the Riemannian volume form of the spatial slices. (This 
equivalence can easily be verified in adapted coordinates.) Thus the integral of pvol over a comoving domain 
remains constant. 



We mention that it is possible to derive the Raychaudhuri equation (4.18) also from the Euler equation 
and the field equation (This is closer to what one does in nonrelativistic fluid dynamics.) The two quite 
different derivations reflect some kind of consistency between field and matter equations. 



As an application of ( 4.18 ) we now show, that there are no static dust solutions in the Newton-Cartan 
theory for A = and that for A > there is just one static solution, which corresponds to the Einstein 
universe. 



By definition a static velocity field u is one with vanishing vorticity, 

(= n af3 d a ® dp) = 0, 

and for which the Lie derivatives of the expansion and the strain vanish: 

l u 9 = 0, L u e J = (e J = e af3 d a ® d ). 



(4.23) 



(4.24) 



Indeed, assume that there is no pressure term in ( 4.19| ), then (4.21) reduces to V u u — 0. Using also the 
staticity conditions in the Raychaudhuri equation ( 4.18p , we find 



4wGp = A - h ap h? a a pG o a(i - -9 2 . 

3 



(4.25) 



This equation has for A = obviously no solution with p > 0. (Note, we have not used the second equation 
of (4.24) to arrive at this conclusion.) 



Consider next the case A > 0. If we write ( 4.25 ) in terms of Galilei coordinates, we obtain 



In such coordinates one has (with the first equation in ( 4.24| )) 

h ij (L u e% = u i d u> !i + ^(i 



(4.26) 



and this vanishes by the second equation of (4.24). Thus the density p satisfies the relation 

A 



P : 



AttG 



of the Einstein universe. 



(4.27) 



These conclusions hold in particular for Newtonian cosmological dust models. It has to be emphasized 
that we have not made any symmetry assumptions (apart from staticity). A very similar argument works 
also in GR ill. 
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5 Newton- Cartan cosmology 



It is very fortunate that the post-recombination universe can be described largely in the Newtonian appro- 
ximation. This brings enormous simplifications in treating the problems of structure formation, in particular 
in the nonlinear regime. Thanks to this circumstance, we can for instance use N-body simulations. 

We consider this as a motivation (beside others) to put Newtonian cosmology on a conceptually firm 
basis. This can readily be achieved in the framework of the geometrical formulation of the Newton-Cartan 
theory that we have described in the previous sections. Again, the analogy to GR is striking. To illustrate 
this, we consider first homogeneous cosmological models and then develop the cosmological perturbation 
theory of Friedmann-Lemaitre models. 



5.1 Homogeneous cosmological models 

In analogy to the discussion of homogeneous cosmological models in GR (for an introduction see we 
consider first the geometrical aspect, without imposing the field equations. Spacetime is then described by 
a Galilei mani fold (M, h, r) with a symmetric Galilei connection T. We introduce adapted coordinates (see 
equations (2.9)). The spatial coordinates {x 1 } parametrize the slices ^ t of constant time on which h induces 
the Riemannian metric h b = hijdx 1 <g> dx 3 . We choose again V = d t . 

Let us assume now that there is a free isometric left action of a 3-dimensional Lie group G on the slices 
Y] t with G on which h defines a time-dependent family of Riemannian metrics. Relative to a left invariant 
basis {9 a } of G this family is of the form h v = h ab (t)9 a ® 9 b . 

Using Vg = V/i = (g = t <E> t = dt ® dt), Cartans structure equations for the connection and the 
Maurer-Cartan equations for the Lie group G, one can then work out the Ricci tensor for all Bianchi types, 
with the result given in |r|. Here, we consider only the Bianchi type I, because the field equations imply 
that the J^t are na *- The metric homogeneity is thus not an additional restriction in the Newton-Cartan 
theory. 

For the choice 9 a = dx a we can compute the Ricci tensor also directly with the help of ( p. 11 ) and set up 
the field equations (3.5). The result is 



Ro 



Roi 



-^(h ij kj),o - -h ij h jk h kl h H + 2h ij K 0jti + K ij Ki 
AnGp - A, 



(5.1) 
(5.2) 



Rij vanishes identically. Equation (5.1) is obviously equivalent to the Raychaudhuri equation (4.18) for 
u = V . The latter reads in adapted coordinates for any velocity field u 

1 



divV„u = V u + t;9 + cr ab a ab - n ab Q ab + AnGp - A. 



(5.3) 



The other field equation (5.2) is equivalent to n iJ j = for u = V, since for any u 



i 



[hacU C , 



hhrU C . 



2^n 



(5.4) 



We give also the expressions for the other kinematical quantities 

1 



e afc 



div(V«M) 



[h aC U C b + h bc U C a + h a b], 
z, 

<a + \h ab h ab , 

U% + 2h ab K 0a ,b + 2h bc K ac U a b + h bc U a b h ac + U a U b ab + U%U b a . 



(5.5) 

(5.6) 
(5.7) 
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Beside this the 2- form k is assumed to be closed (Newtonian connection). 



Matter is assumed to be an ideal fluid with energy momentum tensor (4. IE). In adapted coordinates the 
continuity equation (4.20) and the Euler equation ( [4.21 ) become 



p+(pu i ) ti + -h ij h ij p = 



and 



+ v?u^ + 2h ij n 0j + 2v?{\h ik h kj + h lk n jk ) + -hPpj = 0. 



(5.8) 



(5.9) 



Unlike as in GR, we cannot conclude from our basic equations that the physical quantities like p and p are 
only functions of time. The reason is clear: As already emphasized, the field equations imply that spacetime 
has to be of Bianchi type I. 



Let us specialize the field and matter equations to Newtonian gravity, characterized by condition (3.10). 
We can then introduce Galilei coordinates such that Kij — 0. Relative to these equations ( |5.8| ), (p"ls) and 
3.1|) reduce to 



p + {pu l 



u 1 + u J u l j = -±p ti + g\ 



-g\ = iirGp - A, 



(5.10) 



where g % — ~2n i (as in ( 3.1 if) ) . We thus arrive at the traditional equations for Newtonian gravity, coupled 
to an ideal fluid. 



Let us go back to the Newton-Cartan t heory and assume now that p and p are - in adapted co ordi nates - 
only functions of i. The Euler equation ( 4.21 ) implies then V u u — and the continuity equation (5.8) shows 
that u % i depends only on t. This leads to the following simplification of (5.3) 



1 



a ab a ab - n aD n ab + 4irGp - A = 0. 



Here we have used that 9 is also only a function of t, because (5^) and ( |5.6| ) imply 

p + 6p = 0. 

Specializing again to Newtonian gravity, we can reach stronger conclusions. 



(5.11) 



(5.12) 



Lemma 1 Z/f^+B' is also translational invariant, then there exist for Newtonian gravity Galilei coordinates 
relative to which the spatial components of u — do + u l di are linear functions of the x 3 : 

u l = a){t)x J +b l {t). (5.13) 

Proof: We know that we can introduce Galilei coordinates such that Kij = 0. Since Lg s ,(0" + 0") = 
implies that (O^- + 6y ),fc = 0, equations ( pT5| ) and ( |5.6| ) show that u l j k = 0. This proofs the assertion. 



With a time dependent translation we can pass to Galilei coordinate for which the inhomogeneity in ( 5.13| ) 



disappears. We are now in a situation which has been discussed in classic papers by Heckmann and Schucking 

0, Bi- 



ll 



We consider finally homogeneous and isotropic Friedmann-Lemaitre models in the Newton-Cartan theory. 
We find these with the ansatz 



V = u, 



hij = a 2 (t)Sij, 



= 0. 



(5.14) 



and given by v Vt s . Furthermore, we find 



From the remark connected to equation (2.6) it is clear, that the symmetric Galilei connection is now fixed 



A 

a 



Oab = 0, 



9 = 3- 



and the basic equations (|5.1l| ), ( p.l2|) reduce to 

= 



a 

P + 3-/3, 
a 

AttG 



A 
3"' 



Equation (5.16) implies 



9 = Pa 



oo\ 3 
a) ' 



and when this is used in (5.17) we obtain the Friedmann equation 

■2,, 87rG 2 , A 2 

a +k = —^~pa + —a , 
in which the integration constant k can be choosen to be k = 0, ±1. 



(5.15) 

(5.16) 
(5.17) 

(5.18) 

(5.19) 



In the next section we discuss the perturbation theory of these homogeneous and isotropic solutions in 
the framework of the Newton-Cartan theory. 



5.2 Cosmological perturbation analysis in the Newton-Cartan theory 

We consider cosmological models deviating only by a small amount from a Friedmann-Lemaitre universe, 
which is defined to be the background. Correspondingly we split all geometric and matter variables into 
their background values, indexed by and small deviations 5p, 5p, 5k, etc. 



The Galilei metric (h,r) is kept fixed. This determines the part r, given in (2.5), of the symmetric 
Galilei connection. Because this is just the background connection, we have k 1 -- ^ = 0. The perturbation of 
the connection is entirely described by 5k. We also note that g(u,u) = 1 (t(u) — 1) requires that the four 
velocity field is of the form 



u = do + 5u l di 



(5.20) 



Inserting all this into the field and matter equations leads to a set of perturbation equations for 5p, 5p, 5u l 
and 5k which are still exact. In writing them down, we drop the variational symbol 5 and use the notation 



p = p(°Hl + D) 

In [|l3| the following complete set of perturbation equations is derived: 

i) + [u < (l + D)] li = 0, 



(5.21) 



(5.22) 



it 1 + uV- + 2-u l = -//"/-, - 2// "/,•„, + 2h ij KjiU l 



(5.23) 
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2h ij K 0jli = -K ij Kij + 4irGp w D, 



dn = 0, 



0. 



(5.24) 
(5.25) 
(5.26) 



These agree for k^- = with the usual Newtonian perturbation equations (see, e.g., Jl8[| ) . (In making 
this comparison one has to note, that the peculiar velocity field v l is usually defined by v % = a(t)u l . The 
gravitational field is again given by n oi : g l = —2kqi.) 



Linearization of the perturbation equations gives 

D + u\ = 0, 



u l + 2-u l 
a 



1 



h ij p tj -2h ij K 0j , 



2h ij K 0j ,i = 4irGp {0) D, 



(5.27) 
(5.28) 
(5.29) 



dn = 0, //■/,•,,./ = 0. (5.30) 

Eliminating u l we arrive at the well-known perturbation equation for the density fluctuations: 

D + 2-D = -^h ij p a + AttGp^D. (5.31) 



From (5.22) - (5.26) we can derive in a standard manner (exact) perturbation equations for vorticity and 
shear. One equation agrees with the Raychaudhuri equation for the perturbations: 



AnGp - A = 0. 



For the vorticity one finds E 



(5.32) 



(5.33) 



(We have again dropped the variational symbol on f2, 0, u, p; but p is the total density.) 

We conclude this discussion by writing the exact perturbation equations ( |5.22 ) - ( |5.26] ) in a covariant 



fori 



V V D + div[(l + D){&Z ~ g a xV x V p )u a dp] = 0, 



(5.34) 



V u u — — divph, 
P 



dwV v V = -h^h^K^K-ys + 47rGp (0) L», 



dn = 0, 



ft a,J v p!- Sl /r)^] = o. 



(5.35) 

(5.36) 
(5.37) 
(5.38) 
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6 Concluding remarks 



The Newton-Cartan theory can sometimes provide useful insights for problems in GR. An interesting example 
concerns the cosmic no-hair conjecture, which is not yet settled in sufficient generality within GR. Bauer et 



al 1 19 were, however, able to prove satisfactory theorems in the framework of the Newton-Cartan theory. For 
ideal fluid models they showed that solutions corresponding to nearly homogeneous initial data for a compact 
time slice exist in the case A > for all positive times and that the difference between the inhomogeneous and 
homogeneous solutions tends to zero in a strong sense. Perturbations are thus strongly damped. Presumably 
a corresponding nonlinear stability property holds also in GR, but this appears very difficult to prove. 

The geometrical formulation of the Newton-Cartan theory has also played a useful role in rigorous 
discussions of the Newtonian limit of GR ||. The starting point is the observation by Ehlers that both 
theories fit naturally into a larger frame theory with two metrics h a/3 , g a p related by g QO •ft. <T,3 = — (A = 1 
for GR and A = for the Newton-Cartan theory). 

This frame theory has also played a remarkable role in the work of Heilig |2(J for establishing rigorous 
existence theorems in GR for solutions which describe rotating stars. 

Several more formal aspects have been studied by Kiinzle and collaborators. An example is the genera- 
lization of the Galilei invariant spin-i-wave equation to a curved Newton space-time p| |2l| . 



Finally, without being complete, we mention that Julia and Nicolai 22 have recently obtained the 



Newton-Cartan theory through a dimensional reduction of a Kaluza-Klein theory along a null vector. 

All this demonstrates once more the remarkable continuity in the development of theoretical physics. 
The word "revolution" rarely deserves to be used in this context. To our knowledge, it appears in Einsteins 
writings only once, namely in connection with his hypothesis of the light quantum . He did not regard 
its use to be appropriate in all his work on special and general relativity. 
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